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Stabilization of semilinear PDE’s, and uniform decay under 

discretization 

Emmanuel Trelat* 


Abstract 

These notes are issued from a short course given by the author in a summer school in 
Chambery in June 2015. 

We consider general semilinear PDE’s and we address the following two questions: 

1) How to design an efficient feedback control locally stabilizing the equation asymptotically 
to 0? 

2) How to construct such a stabilizing feedback from approximation schemes? 

To address these issues, we distinguish between parabolic and hyperbolic semilinear PDE’s. 
By parabolic, we mean that the linear operator underlying the system generates an analytic 
semi-group. By hyperbolic, we mean that this operator is skew-adjoint. 

We first recall general results allowing one to consider the nonlinear term as a perturbation 
that can be absorbed when one is able to construct a Lyapunov function for the linear part. 
We recall in particular some known results borrowed from the Riccati theory. 

However, since the numerical implementation of Riccati operators is computationally de¬ 
manding, we focus then on the question of being able to design “simple” feedbacks. For 
parabolic equations, we describe a method consisting of designing a stabilizing feedback, based 
on a small finite-dimensional (spectral) approximation of the whole system. For hyperbolic 
equations, we focus on simple linear or nonlinear feedbacks and we investigate the question of 
obtaining sharp decay results. 

When considering discretization schemes, the decay obtained in the continuous model 
cannot in general be preserved for the discrete model, and we address the question of adding 
appropriate viscosity terms in the numerical scheme, in order to recover a uniform decay. We 
consider space, time and then full discretizations and we report in particular on the most 
recent results obtained in the literature. 

Finally, we describe several open problems and issues. 
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Universitaire de France, F-75005, Paris, France (enraianuel.trelatOupnic.fr). 
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1 Introduction and general results 

1.1 General setting 

Let X and U be Hilbert spaces. We consider the semilinear control system 

y{t) = Ay{t)+F(y(t)) + Bu(t), (1) 

where A : D(A ) —► X is an operator on X , generating a Co semi-group, B £ L(U,D(A*)'), and 
F : X —> X is a (nonlinear) mapping of class C 1 , assumed to be Lipschitz on the bounded sets 
of X , with F(0) = 0. We refer to [10, 18, 43, 44] for wcll-posedness of such systems (existence, 
uniqueness, appropriate functional spaces, etc). 

We focus on the following two questions: 

1. How to design an efficient feedback control u = Ky , with K £ L(X,U), locally stabilizing 
(1) asymptotically to 0? 

2. How to construct such a stabilizing feedback from approximation schemes? 

Moreover, we want the feedback to be as simple as possible, in order to promote a simple imple¬ 
mentation. Given the fact that the decay obtained in the continuous setting may not be preserved 
under discretization, we are also interested in the way one should design a numerical scheme, in 
order to get a uniform decay for the solutions of the approximate system, i.e., in order to guaran¬ 
tee uniform properties with respect to the discretization parameters Ax and/or At. This is the 
objective of these notes, to address those issues. 

Concerning the first point, note that, in general, stabilization cannot be global because there 
may exist other steady-states than 0, i.e., y £ X such that Ay + F(y) = 0. This is why we speak 
of local stabilization. Without loss of generality we focus on the steady-state 0 (otherwise, just 
design a feedback of the kind u = K(y — y)). 

We are first going to recall well known results on how to obtain local stabilization results, 
first in finite dimension, and then in infinite dimension. As a preliminary remark, we note that, 
replacing if necessary A with A + dF( 0), we can always assume, without loss of generality, that 
dF( 0) = 0, and thus, 

\\ F (y)\\ = o(IMI) near y = 0. 

Then, a first possibility to stabilize (1) locally at 0 is to consider the nonlinear term F(y) as a 
perturbation, that we are going to absorb with a linear feedback. Let us now recall standard results 
and methods. 

1.2 In finite dimension 

In this subsection, we assume that X = 1R”. In that case, A is a square matrix of size n, and 
B is a matrix of size n x m. Then, as it is well known, stabilization is doable under the Kalman 
condition 

rank(H, AB ,..., A n ~ l B) = n, 

and there are several ways to do it (see, e.g., [42] for a reference on what follows). 
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First possible way: by pole-shifting. According to the pole-shifting theorem , there exists K 
(matrix of size m x n) such that A + BK is Hurwitz, that is, all its eigenvalues have negative real 
part. Besides, according to the Lyapunov lemma , there exists a positive definite symmetric matrix 
P of size n such that 

P(A + BK) + (A + BK) t P = -I n . 

It easily follows that the function V defined on 1R” by 

V{y) = y T Py 

is a Lyapunov function for the closed-loop system y(t) = (A + BK)y(t). 

Now, for the semilinear system (1) in closed-loop with u[t) = Ky[t), we have 

J t v (y(t)) = -hm 2 + y{t) T PF{y(t)) < -C4y{t)f < - C 2 V(y(t )), 

under an a priori assumption on y(t), for some positive constants C\ and C 2 , whence the local 
asymptotic stability. 

Second possible way: by Riccati procedure. The Riccati procedure consists of computing 
the unique negative definite solution of the algebraic Riccati equation 

A t E + EA + EBB T E = I n . 


Then, we set u(t) = B T Ey(t). 

Note that this is the control minimizing / 0 +OO (||2/(i)|| 2 + ||M(t)|| 2 )c?t for the control system y(t) = 
Ay{t) + Bu(t). This is one of the best well known results of the Linear Quadratic theory in optimal 
control. 

Then the function 

V{y) = —y*Ey 

is a Lyapunov function, as before. 

Now, for the semilinear system (1) in closed-loop with u(t) = B T Ey(t), we have 
jV(y{t)) = -y(t) r (I n + EBB T E)y{t) - y(t) T EF(y(t)) < -C^yWf < -C 2 V(y(t)), 
under an a priori assumption on y{t ), and we easily infer the local asymptotic stability property. 

1.3 In infinite dimension 

1.3.1 Several reminders 

When the Hilbert space X is infinite-dimensional, several difficulties occur with respect to the 
finite-dimensional setting. To explain them, we consider the uncontrolled linear system 

y{t) = Ay(t), (2) 

with A : D(A) —► X generating a Co semi-group S(t). 

The first difficulty is that none of the following three properties are equivalent: 

• S(t) is exponentially stable, i.e., there exist C > 0 and S > 0 such that ||<S'(t)|| < Ce~ St , for 
every t ^ 0. 

• the spectral abscissa is negative, i.e., sup{Re(A) | A € a (A)} < 0. 
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• All solutions of (2) converge to 0 in X , i.e., S(t)y° —> 0, for every y° £ X. 

£->-+oo 

For example, if we consider the linear wave equation with local damping 


ytt - A y + Xu,y = 0, 

in some domain fl of JR/ 1 , with Dirichlet boundary conditions, and with w an open subset of G, then 
it is always true that any solution ( y,y t ) converges to (0,0) in Hq(Q) x L 2 (G) (see [16]). Besides, 
it is known that we have exponential stability if and only if w satisfies the Geometric Control 
Condition (GCC). This condition says, roughly, that any generalized ray of geometric optics must 
meet ui in finite time. Hence, if for instance fl is a square, and u is a small ball in f2, then GCC 
does not hold and hence the exponential stability fails, whereas convergence of solutions to the 
equilibrium is valid. 

In general, we always have 

sup{Re(A) | A £ a(A)} ^ inf{/z £ H | 3C > 0, ||S(i)|| < Cef* Vt > 0}, 

in other words the spectral abscissa is always less than or equal to the best exponential decay rate. 
The inequality may be strict, and the equality is referred to as “spectral growth condition”. 

Let us go ahead by recalling the following known results: 

• Datko theorem: S(t) is exponentially stable if and only if, for every y° £ X , S(t)y° converges 
exponentially to 0. 

• Arendt-Batty theorem: If there exists M > 0 such that ||5(£)|| ^ M for every t ^ 0, and if 

C p(A) (where p(A) is the resolvent set of A), then S(t)y° —> 0 for every y° £ X. 

£—>•+00 

• Huang-Pruss theorem: Assume that there exists M > 0 such that ||S(t)|| ^ M for every t ^ 0. 
Then S(t) is exponentially stable if and only if dR C p(A) and sup^g^ ||(*/3id —A) _1 || < +oo. 

Finally, we recall that: 

• Exactly null controllable implies exponentially stabilizable, meaning that there exists K £ 
L(X, U) such that A + BK generates an exponentially stable Co semi-group. 

• Approximately controllable does not imply exponentially stabilizable. 

For all reminders done here, we refer to [15, 18, 31, 32, 45]. 

1.3.2 Stabilization 

Let us now consider the linear control system 

y = Ay + Bu 

and let us first assume that B £ L(U,X), that is, the control operator B is bounded. We assume 
that the pair (A, B) is exponentially stabilizable. 


5 


Riccati procedure. As before, the Riccati procedure consists of finding the unique negative 
definite solution E £ L(X) of the algebraic Riccati equation 

A*E + EA + EBB*E = id, 

in the sense of ((2 EA + EBB*E — id )y,y) = 0, for every y £ D{A), and then of setting u(t) = 
B*Ey(t). Note that this is the control minimizing / 0 + °°(||2/(t)|| 2 + ||it(i)[| 2 ) dt for the control system 
y = Ay + Bu. 

Then, as before, the function V(y) = —{y, Ey) is a Lyapunov function for the system in closed- 
loop y = {A + BK)y. 

Now, for the semilinear system (1) in closed-loop with u(t) = B*Ey(t), we have 
jV{y{i)) = (id + EBB*E)y{t)) - ( y(t ), EF(y(t))) < -C? x ||j,(t)|| 2 < -C 2 V(y(t)) 

and we infer local asymptotic stability. 

For B £ L(U,D(A*Y ) unbounded, things are more complicated. Roughly, the theory is com¬ 
plete in the parabolic case (i.e., when A generates an analytic semi-group), but is incomplete in 
the hyperbolic case (see [29, 30] for details). 

Rapid stabilization. An alternative method exists in the case where A generates a group S(t), 
and B £ L(U, D(A*)') an admissible control operator (see [44] for the notion of admissibility). An 
example covered by this setting is the wave equation with Dirichlet control. 

The strategy developed in [27] consists of setting 

rT+ 1/2A 

u = -B*Cy 1 y with C x = fx{t)S(—t)BB*S(—t)* dt, 

J 0 

with A > 0 arbitrary, f\{t) = e~ 2Xt iff £ [0, T] and fx(t) = 2\e~ 2XT (T+1/2X—t) if t £ [T, T+1/2A]. 
Besides, the function 

V(y) = (y,Cy 1 y) 

is a Lyapunov function (as noticed in [11]). Actually, this feedback yields exponential stability with 
rate —A for the closed-loop system y = (A— BB*C^ 1 )y, whence the wording “rapid stabilization” 
since A > 0 is arbitrary. 

Then, thanks to that Lyapunov function, the above rapid stabilization procedure applies as 
well to the semilinear control system (1), yielding a local stabilization result (with any exponential 
rate). 

1.4 Existing results for discretizations 

We recall hereafter existing convergence results for space semi-discretizations of the Riccati proce¬ 
dure. We denote by En the approximate Riccati solution, expecting that Ejy —> E as Af —► +oo. 

One can find in [7, 8, 22, 26, 32] a general result showing convergence of the approximations 
Ejsr of the Riccati operator E , under assumptions of uniform exponential stabilizability, and of 
uniform boundedness of En". 

ll‘5 , A A r+BNB^B N ( i )ll ^ Ce St i PwKM, 

for every t ^ 0 and every N £ IN*, with uniform positive constants C, 5 and M. 
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In [7, 29, 32], the convergence of Em to E is proved in the general parabolic case, for un¬ 
bounded control operators, that is, when A : D(A) —> X generates an analytic semi-group, 
B £ L(U, D(A*)'), and (A,B) is exponentially stabilizable. 

The situation is therefore definitive in the parabolic setting. In contrast, if the semi-group 
S(t) is not analytic, the theory is not complete. Uniform exponential stability is proved under 
uniform Huang-Priiss conditions in [32], More precisely, it is proved that, given a sequence (£„(•)) 
of Cq semi-groups on X n . of generators A n . (S n (-)) is uniformly exponentially stable if and only if 
ilR C p(A n ) for every n £ IN and 

sup || (*/3id - A n y 1 \\ < + 00 . 

/3eR, nelN 

This result is used e.g. in [37] to prove uniform stability of second-order equations with (bounded) 
damping and with viscosity term, under uniform gap condition on the eigenvalues. 

This result also allows to obtain convergence of the Riccati operators, for second-order systems 

y + Ay = Bu , 

with A : D(A ) — > X positive selfadjoint, with compact inverse, and B £ L(U,X) (bounded control 
operator). 

But the approximation theory for general LQR problems remains incomplete in the general 
hyperbolic case with unbounded control operators, for instance it is not done for wave equations 
with Dirichlet boundary control. 

1.5 Conclusion 

Concerning implementation issues, solving Riccati equations (or computing Gramians, in the case 
of rapid stabilization) in large dimension is computationally demanding. In what follows, we would 
like to find other ways to proceed. 

Our objective is therefore to design simple feedbacks with efficient approximation procedures. 
In the sequel, we are going to investigate two situations: 

• Parabolic case ( A generates an analytic semi-group): we are going to show how to design 
feedbacks based on a small number of spectral modes. 

• Hyperbolic case, i.e., A = — A* in (1): we are going to consider two “’’simple” feedbacks: 

— Linear feedback u = — B*y : in that case, if F = 0 then 5 ^||y(i)|| 2 = — \\B*y(t)\\ 2 . We 
will investigate the question of how to ensure uniform exponential decay for approxima¬ 
tions. 

— Nonlinear feedback u = B*G(y ): we will as well investigate the question of how to 
ensure uniform (sharp) decay for approximations. 


2 Parabolic PDE’s 

In this section, we assume that the operator A in (1) generates an analytic semi-group. Our 
objective is to design stabilizing feedbacks based on a small number of spectral modes. 

To simplify the exposition, we consider a ID semilinear heat equation, and we will comment 
further on extensions. 

Let L > 0 be fixed and let / : HR —>• 1R be a function of class C 2 such that /(0) = 0. Following 
[12], We consider the ID semilinear heat equation 

yt = yxx + f{y), y{t,0) = 0, y{t,L) = u(t), (3) 
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where the state is y(t, •) : [0, L\ —>• IR and the control is u(t) G H. 

We want to design a feedback control locally stabilizing (3) asymptotically to 0. Note that this 
cannot be global, because we can have other steady-states (a steady-state is a function y G C 2 ( 0, L) 
such that y"(x) + f(y(x)) = 0 on (0, L ) and y(0) = 0). By the way, here, without loss of generality 
we consider the steady-state 0. 

Let us first note that, for every T > 0, (3) is well posed in the Banach space 

Y t = L 2 { 0, T; H 2 { 0, L)) fl H\ 0, T; L 2 (0, L)), 

which is continuously embedded in L°°((0,T) x (0,-L)). 1 

First of all, in order to end up with a Dirichlet problem, we set 

z(t,x) = y(t, x) - ^ u(t ). 

Assuming (for the moment) that u is differentiable, we set v(t) = u'(t ), and we consider in the 
sequel v as a control. We also assume that it(0) = 0. Then we have 

Zt = Zxx + f'(0)z + - ^-v + r(t,x), z(t,0) = z(t,L) = 0, (4) 

with z(0,x) = 2 /( 0 , cc) and 


r(t,x) = (z(t,x) + ^u{t)^ J (1 — s)f" ^sz(s, x) + s—u(s)j ds. 

Note that, given B > 0 arbitrary, there exist positive constants C\ and C 2 such that, if |«(t)| < B 
and || z(t, ■) ||l°°(o,l) < B, then 

IKV)IIl~( 0,L) < + ||z(t,-)llL»(0,£)) < C 2 (u{t ) 2 + |k(t,-)llffi( 0 ,L))- 


In the sequel, r(t,x) will be considered as a remainder. 

We define the operator A = A + /'(0)id on D{A) = H 2 { 0, L) fl Hq (0, L), so that (4) is written 
as 

u = v , Zt = Az + au + bv + r, z(t, 0) = z(t, L) = 0, (5) 

with a(x) = f f'(0) and b(x) = -f. 

Since A is selfadjoint and has a compact resolvent, there exists a Hilbert basis 1 of 

L 2 (0,L), consisting of eigenfunctions e.j G Hq(0,L) fl C 2 ([0,L]) of A , associated with eigenvalues 
such that —00 < • • • < \ n < ■ ■ ■ < Ai and \ n —> —00 as n —> + 00 . 

Any solution z(t, •) G H 2 ( 0, L) fl Hq (0, L) of (4), as long as it is well defined, can be expanded 
as a series 

OO 

z (tr) 

j =1 


1 Indeed, considering v G L 2 (0,T; H 2 (0, L)) with vt G H 1 (0, T; L 2 (0, L)), writing v = Cjke’’i t e’' kx , we have 


j,k 


^I 


+ j 2 + k 4 


1/2 


TV + J 2 + k 4 )\cjk\' 


1/2 


and these series converge, whence the embedding, allowing to give a sense to f(y). 

Now, if yi and y 2 are solutions of (3) on [0, T], then yi = y 2 - Indeed, v = y\ — y 2 is solution of vt = v X x fan, 
v(t, 0) = v(t, L) = 0, u(0,ai) = 0, with a(t,x) = g(yi(t i x) i y 2 (t 1 x)) where g is a function of class C 1 . We infer that 
v = 0. 
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(converging in Hq(0,L)), and then we have, for every j ^ 1, 


Zj(t) = A jZj{t) + a,ju(t) + bj V (t) + rj{t ), 


with 


m 


r*L />L /*L 

/ xej(x)dx, bj=—— / xej(x)dx, fj(i) = / r(t,x)ej(x)dx. 

Jo L Jo Jo 


Setting, for every n £ IN* 


X n (t) = 


we have, then, 


/ u(t) \ 


(0 

o • 



f l ) 


f 0 \ 

Zl(t) 


0.1 

Ai 

• 0 


h 

, R n {t) = 

n(t) 

> A n = 




, Bn = 


\Zn{t)j 


\P*n 

0 • 

• aJ 


\b n j 


Vn{t)j 


Xn(t) — A n X n (t) + B n v(t ) + R n (t). 
Lemma 1. The pair ( A ni B n ) satisfies the Kalman condition. 
Proof. We compute 


det (B n , A n B n ,..., A™B n ) = JJ(aj +A jbj) VdM(Ai,..., A„), 

j =i 


( 6 ) 


where VdM(Ai,..., A ra ) is a Van der Monde determinant, and thus is never equal to zero since 
the A i, i = 1... n, are pairwise distinct. On the other part, using the fact that each ej is an 
eigenfunction of A and belongs to Hq (0, L), we compute 

1 f L i f L 

a j + A jbj = — x(f'( 0) — Aj)ej(x) dx = —— / xe'Ux ) dx = — e'AL), 

L Jo L Jo 

and this quantity is never equal to zero since Oj(L) = 0 and ej is a nontrivial solution of a linear 
second-order scalar differential equation. Therefore the determinant ( 6 ) is never equal to zero. □ 

By the pole-shifting theorem, there exists K n = (fco,..., k n ) such that the matrix A n + B n K n 
has —1 as an eigenvalue of multiplicity n + 1. Moreover, by the Lyapunov lemma, there exists a 
symmetric positive definite matrix P n of size n + 1 such that 

Pn ( A n + B n K n ) + (A n + B n K n ) P n = —I n +\. 

Therefore, the function defined by V n (X) = X T P n X for any X £ H ™" 1 " 1 is a Lyapunov function 
for the closed-loop system X n (t) = (A n + B n K n )X n (t). 

Let 7 > 0 and n £ IN* to be chosen later. For every u £ 1R and every z £ H 2 { 0, L) (~l (0, L), 

we set 

1 1 °° 

V(u, z) = 7 XlP n X n - -(z, Az) L 2 { 0iL) = 7 XjP n X n - - ^ A^|, (7) 

j=i 

where X n £ lR n+1 is defined by X n = (u, z i,.. •, z n ) T and Zj = (z(-), et{-)) l 2 (o,l) for every j. 
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Using that \ n —> —oo as n —V +oo, it is clear that, choosing 7 > 0 and n £ IN* large enough, 
we have V (it, z) > 0 for all (u, z) £ R x (H 2 ( 0, L) (~l Hq(0, L)) \ {(0,0)}. More precisely, there exist 
positive constants C3, C4, C5 and Cq such that 

C 3 (u +11^11^1(0^)) ^ V{ u , z ) ^ C 4 (u + INIffifO,!,)) ) 

V(u,z) < C 6 (\\X n \\l + \\Az\\h l0 ,L)), lC 6 \\X n \\l < V(u,z), 

for all (u, z) £ R x (ff 2 (0, L) D Hq(0, L)). Here, || H 2 designates the Euclidean norm of 1R" +1 . 

Our objective is now to prove that V is a Lyapunov function for the system (5) in closed-loop 
with the control v = K n X n . 

In what follows, we thus take v = K n X n and u defined by u = v and u(0) = 0. We compute 

^V(u(t),z(t)) = - 7 \\X n (t)\\l - \\Az(t,-)f L 2 - (Az(t r ),a(-)) L 2 U (t) 

- ( Az(t , •), b(-)) L 2 K n X n (t) - (Az{t , •), r(t, -)} L 2 

+ 7 ( R n (t) T P n X n (t ) + X n (t) T P n R n (t )) . ( 8 ) 

Let us estimate the terms at the right-hand side of ( 8 ). Under the a priori estimates \u(t)\ ^ B 
and \\z(t, -)||.l«>(o,£) ^ B, there exist positive constants C 7 , Cs and Cg such that 

I (Az, a) L 2u\ + \{Az, b) L 2 K n X n \ < -||H2||| 2 + C'rll-X'nlll, 

\(Az,r) L 2\ < l\\Az \\ 2 L2 + c 8 v 2 , Halloo < §V, 

|7 {RlP n X n +XjP n R n ) | < J^—^V*/ 2 . 

03VU 

We infer that, if 7 > 0 is large enough, then there exist positive constants C 10 and C\\ such that 
^V < —CiqV + CnU 3 / 2 . We easily conclude the local asymptotic stability of the system (5) in 
closed-loop with the control v = K n X n . 

Remark 1 . Of course, the above local asymptotic stability may be achieved with other procedures, 
for instance, by using the Riccati theory. However, the procedure developed here is much more 
efficient because it consists of stabilizing a finite-dimensional part of the system, mainly, the part 
that is not naturally stable. We refer to [12] for examples and for more details. Actually, we have 
proved in that reference that, thanks to such a strategy, we can pass from any steady-state to any 
other one, provided that the two steady-states belong to a same connected component of the set 
of steady-states: this is a partially global exact controllability result. 

The main idea used above is the following fact, already used in the remarkable early paper 
[38]. Considering the linearized system with no control, we have an infinite-dimensional linear 
system that can be split, through a spectral decomposition, in two parts: the first part is finite¬ 
dimensional, and consists of all spectral modes that are unstable (meaning that the corresponding 
eigenvalues have nonnegative real part); the second part is infinite-dimensional, and consists of all 
spectral modes that are asymptotically stable (meaning that the corresponding eigenvalues have 
negative real part). The idea used here then consists of focusing on the finite-dimensional unstable 
part of the system, and to design a feedback control in order to stabilize that part. Then, we plug 
this control in the infinite-dimensional system, and we have to check that this feedback indeed 
stabilizes the whole system (in the sense that it does not destabilize the other infinite-dimensional 
part). This is the role of the Lyapunov function V defined by (7). 
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The extension to general systems (1) is quite immediate, at least in the parabolic setting under 
appropriate spectral assumptions (see [39] for Couette flows and [14] for Navier-Stokes equations). 

But it is interesting to note that it does not work only for parabolic equations: this idea has 
been as well used in [13] for the ID semilinear equation 

ytt = Vxx + f{y), y(t, o) = o, y x (t , l) = u(t), 

with the same assumptions on / as before. We first note that, if f(y) = cy is linear (with 
c € L°°(0,L)), then, setting u(t) = — ayt(t,L ) with a > 0 yields an exponentially decrease of the 
energy (y t (t,x) 2 + y x (t,x) 2 ) dt, and moreover, the eigenvalues of the corresponding operator 
have a real part tending to —oo as a tends to 1. Therefore, in the general case, if a is sufficiently 
close to 1 then at most a finite number of eigenvalues may have a nonnegative real part. Using 
a Riesz spectral expansion, the same kind of method as the one developed above can therefore 
be applied, and yields a feedback based on a finite number of modes, that stabilizes locally the 
semilinear wave equation, asymptotically to equilibrium. 


3 Hyperbolic PDE’s 

In this section, we assume that the operator A in (1) is skew-adjoint, that is, 

A*=-A, D{A*) = D{A). 

Let us start with a simple remark. If F = 0 (linear case), then, choosing the very simple linear 
feedback u = —B*y and setting V(y) = IgjHj/Hx, we have 

lv(y(t)) = ~\\B*y(t)\\ 2 x ^0, 

and then we expect that, under reasonable assumptions, we have exponential asymptotic stability 
(and this will be the case under observability assumptions, as we are going to see). 

Now, if we choose a nonlinear feedback u = B*G(y), we ask the same question: what are 
sufficient conditions ensuring asymptotic stability, and if so, with which sharp decay? 

Besides, we will investigate the following important question: how to ensure uniform properties 
when discretizing? 

3.1 The continuous setting 

3.1.1 Linear case 

In this subsection, we assume that F = 0 (linear case), and we assume that B is bounded. Taking 
the linear feedback u = — B*y as said above, we have the closed-loop system y = Ay — BB*y. 
For convenience, in what follows we rather write this equation in the form (more standard in the 
literature) 

y{t) + Ay(t) + By(t) = 0, (9) 

where A is a densely defined skew-adjoint operator on X and B is a bounded nonnegative selfadjoint 
operator on X (we have just replaced A with — A and BB* with B). 

We start hereafter with the question of the exponential stability of solutions of (9). 
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Equivalence between observability and exponential stability. The following result is a 
generalization of the main result of [23] . 

Theorem 1. Let X be a Hilbert space, let A : D(A) —» X be a densely defined skew-adjoint 
operator, let B be a bounded selfadjoint nonnegative operator on X. We have equivalence of: 

1. There exist T > 0 and C > 0 such that every solution of the conservative equation 

4>(t) + A(j>{t) = 0 (10) 


satisfies the observability inequality 


wmfx^c f 0 \\B^m\\dt. 

Jo 


2. There exist C 1 > 0 and 5 > 0 such that every solution of the damped equation 


y(t ) + Ay(t ) + By(t) = 0 


( 11 ) 


satisfies 

E y {t) < CiE y (0)e- ft , 

where E y (t) = \\\y(t)\\ 2 x . 

Proof. Let us first prove that the first property implies the second one: we want to prove that 
every solution of ( 11 ) satisfies 


E v (t) = \hmx < E y (0)e~ st = i|| 2 /( 0 )||ie- 5t . 

Consider <\> solution of (10) with </>(0) = y( 0). Setting 0 = y — we have 

6 + A6 + By = 0, 0(0) = 0. 

Then, taking the scalar product with 0, since A is skew-adjoint, we get (0 + By,6)x = 0. But, 
setting E e (t) = \\\8{t)\\ 2 x , we have Eg = — ( By,6)x■ Then, integrating a first time over [0,i], and 
then a second time over [0, T], since Eg(0) = 0, we get 

[ Eg(t)dt = - I [ {By(s),0(s)) x dsdt = - [ {T - t){B 1 / 2 y{t),B 1 / 2 6(t)) x dt, 

Jo Jo Jo Jo 

where we have used the Fubini theorem. Hence, thanks to the Young inequality ab ^ ^a 2 + ^fe 2 
with a = 2 , we infer that 


\ £ wmfx dt < T\\B l / 2 \\ £ \\B l / 2 y(f)\\x\\m\\x dt 

< T 2 ||H 1/2 || 2 £ \\B l / 2 y{t)\\ 2 x dt+\£ \m\\x dt, 

and therefore, 

£ \m\\x dt < 4T 2 ||H 1 / 2 |||- £ \\B^ 2 y{t)\\x dt. 

Jo Jo 
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Now, since (f> = y — 9, it follows that 



\ Bl/2 Ht)\\xdt^2 f \\B 1/2 y(t)\\ 2 x dt + 2 f \\B 1/2 d(t)\\ 2 x dt 
Jo Jo 


^ (2 + 8T 2 \\B 1/2 \\ A 


IB^ymldt. 


Finally, since 

Ey(0) = E^0)= 1 -\\m\\x < T £ WB^m&dt 


it follows that 


Ey( 0) < CO + 4T 2 || j B 1 / 2 || 4 ) [ T \\B^y{t)f x dt. 

Jo 

Besides, one has E' y (t ) = — \\B 1 / 2 y{t )\\and then / 0 T \\B 1 / 2 y(t)\\ 2 x dt = E y { 0) — E y {T). Therefore 
Ey{ 0) < C(1 + 4T 2 ||B 1 / 2 || 4 )(£ , y (0) - Ey(T)) = <7^(0) - E y (T )) 


and hence 

E v (T)^^-J E y(0) = C 2 Ey(0), 

with Ci < 1. 

Actually this can be done on every interval [ kT , (fc-fT)T], and it yields E y ((k+1)T) ^ C 2 E y (kT) 
for every k € IN, and hence E y (kT) ^ £ , y (0)C'f. 

For every t £ [kT, (k + 1 )T), noting that k = [y] > y — 1, and that In ^ > 0, it follows that 

Cj = exp(fclnC 2 ) = exp(—fc In —!p) ^ J- exp(— 

02 02 

and hence E y (t) < E y (kT) < SE y ( 0) exp(— St) for some 6 > 0. 

Let us now prove the converse: assume the exponential decrease of solutions of (11), and let 
us prove the observability property for solutions of (10). 

From the exponential decrease inequality, one has 

[ T \\B 1/2 y(t)\\ 2 x dt = Ey{ 0) - E y (T) > (1 - C ie - ST )E y ( 0) = C 2 Ey( 0), (12) 

Jo 

and for T > 0 large enough there holds C 2 = 1 — Cie -<5T > 0. 

Then we make the same proof as before, starting from (10), that we write in the form 

</> + A(f> + B(f> = B(j>, 


and considering the solution of (11) with y( 0) = <(>(0). Setting 9 = cj) — y, we have 

6 + Ad + B9 = Bcj), 6>(0) = 0. 

Taking the scalar product with 9, since A is skew-adjoint, we get (9 + B9, 9) x = {B<f>,9) x , and 
therefore 

Eg + (B9,e) x = (Bcj>,e)x. 
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Since ( B6,6)x = ||i? 1//2 6>||^ > 0, it follows that Eg < {B<j>,9)x- As before we apply f^fg and 
hence, since Eg( 0) = 0, 


E e (t)dt < [ f {B(f){s),e(s))xdsdt= [ (T - f)(S 1/2 <(>(f), B 1 d 2 6(t)) x dt. 
Jo Jo Jo 


Thanks to the Young inequality, we get, exactly as before, 


imWxdt^TWB^W / \\B^m\x\\e{t)\\xdt 


< r 2 ii5 1/2 iix £ \\B i/2 m\x dt +\[ M^Wx 


and finally, 

f T ||0(t)|| 2 Y dt < AT*\\B^fx [ T \\B 1/2 m\ 2 x dt. 

JO Jo 

Now, since y = (j) — 0, it follows that 


[ T \\B 1/2 y(t)\\x dt s? 2 f \\B^m\ 2 xdt +2 P WB^O^Wx dt 
Jo Jo Jo 

< (2 + 8T 2 ||R 1/2 || 4 ) [ II B l / 2 cf>{t)\\ 2 x dt 

Jo 

Now, using (12) and noting that -£^(0) = E^O), we infer that 

C 2 E^{ 0) < (2 + 8T 2 ||B 1//2 || 4 ) [ WB^^Ofxdt. 

Jo 


This is the desired observability inequality. 


□ 


Remark 2. This result says that the observability property for the linear conservative equation 
(10) is equivalent to the exponential stability property for the linear damped equation (11). This 
result has been written in [23] for second-order equations, but the proof works exactly in the same 
way for more general first-order systems, as shown here. More precisely, the statement done in [23] 
for second-order equations looks as follows: 

We have equivalence of: 

1. There exist T > 0 and C > 0 such that every solution of 

<p(t) + A<j)(t) = 0 (conservative equation) 

satisfies 

\\A^m\\x + ii ^( o ) ii 2 y < c f To wB^mwdt, 

Jo 

2. There exist C± > 0 and S > 0 such that every solution of 

y(t) + Ay{t) + By(t) = 0 (damped equation) 

satisfies 

E y (t ) < CiT; y (0)e- 5t , 

where 

E y (t) = l(\\A i/ Mm 2 x + \\m\ 2 x)- 
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Remark 3. A second remark is that the proof uses in a crucial way the fact that the operator B is 
bounded. We refer to [5] for a generalization for unbounded operators with degree of unboudedness 
^ 1/2 (i.e., B £ L(U, D(A 1 ' 2 )')), and only for second-order equations, with a proof using Laplace 
transforms, and under a condition on the unboundedness of B that is not easy to check (related 
to “hidden regularity” results), namely, 

V/3 > 0 sup \\B*X(X 2 I + A)~ 1 B\\ LI jj) < +oo. 

Re(A)=/3 


For instance this works for waves with a nonlocal operator B corresponding to a Dirichlet condition, 
in the state space L 2 x H~ x , but not for the usual Neumann one, in the state space H 1 x L 2 (except 
in ID). 


3.1.2 Semilinear case 

In the case with a nonlinear feedback, still in order to be in agreement with standard notations 
used in the existing literature, we rather write the equation in the form u + Au + F(u) = 0, where 
u now designates the solution (and not the control). 

Therefore, from now on and throughout the rest of the paper, we consider the differential system 

u r (t) + Au(t) + B F = 0, (13) 

with A : D{A) — > X a densely defined skew-adjoint operator, B : X — > X a nontrivial bounded 
selfadjoint nonnegative operator, and F:I-)Ia (nonlinear) mapping assumed to be Lipschitz 
continuous on bounded subsets of X. This is the framework and notations adopted in [4]. 

If F = 0 then the system (13) is purely conservative, and ||w(f)||x = ||u(0)||x for every t > 0. 
If F ^ 0 then the system (13) is expected to be dissipative if the nonlinearity F has “the good 
sign”. Along any solution of (13) (while it is well defined), the derivative with respect to time of 
the energy E u (t) = |||u(t)||x is 

E' u (t) = ~{u(t),BF(u{t)))x = —(B 1 / 2 u(t)iB 1 / 2 F{u{t)))x ■ 


In the sequel, we will make appropriate assumptions on B and on F ensuring that E' u (t ) ^ 
0. It is then expected that the solutions are globally well defined and that their energy decays 
asymptotically to 0 as t —> +oo. 

We make the following assumptions. 


• For every u £ X 


(u,BF(u)) x > 0 . 


This assumption implies that E' u (t) ^ 0. 


The spectral theorem applied to the bounded nonnegative selfadjoint operator B implies that 
B is unitarily equivalent to a multiplication: there exist a probability space (fl,p.), a real-valued 
bounded nonnegative measurable function b defined on X , and an isometry U from L 2 (fl,/r) into 
X , such that U~ 1 BUf = bf for every / £ i 2 (fl,/r). 

Now, we define the (nonlinear) mapping p : L 2 (f2,/z) —>- L 2 (f2, p) by 

p(f) = U~ 1 F(U /). 


We make the following assumptions on p: 

• p{ 0) = 0 and fp{f) ^ 0 for every / £ L 2 (fl, p). 
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• There exist Ci > 0 and C 2 > 0 such that, for every / £ L°°(n, p), 

cig(\f(x)\) ^ \p(f)(x)\ s? c 2 3 _ 1 (|/(a:)|) for almost every x £ fl such that |/(x)| < 1, 

c i \f( x )\ ^ \p(f)i x )\ ^ c 2 |/(cc)| for almost every x £ 12 such that |/(x)| ^ 1 , 

where g is an increasing odd function of class C 1 such that 5 ( 0 ) = g'(0) = 0, sg'(s) 2 * /g(s) —> 0 
as s —> 0, and such that the function H defined by H(s) = y/sg(y/s), for every s £ [0,1], is 
strictly convex on [0, s(j] for some so £ (0,1]. 

This assumption is issued from [2] where the optimal weight method has been developed. 
Examples of such functions g are given by 

g(s) = s/ln p (l/s), s p , e -1//s , s p ln 9 (l/s), 

We define the function H on H by H(s) = H(s) for every s £ [0,Sq] and by H(s) = +00 

otherwise. We define the function L on [0, +00) by L(0) = 0 and, for r > 0, by 

L(r) = ^ ^ = - sup (rs — H(s)] , 
r r sgIR V / 


where H* is the convex conjugate of H. By construction, the function L : [0, + 00 ) [0, Sq) is 

continuous and increasing. 

We define A h ■ (0, Sq] —> (0, + 00 ) by A h(s) = H(s)/sH'(s), and we set 


Vs > 1 /H'(sl) 


V’(s) 


1 

HW) 



« 2 (1-A dv - 


The function i/’ : [\/H'{sq), + 00 ) —>■ [0,+oo) is continuous and increasing. 

Hereafter, we use the notations < and — in the estimates, with the following meaning. Let S 
be a set, and let F and G be nonnegative functions defined on 1R x £1 x S. The notation F < G 
(equivalently, G > F) means that there exists a constant C > 0, only depending on the function g 
or on the mapping p, such that F(t, x, A) ^ CG(t, x, A) for all (t,x, A) GlRxfix S. The notation 
Fi ~ F -2 means that F\ < F 2 and F\ > F 2 . 

In the sequel, we choose S = X, or equivalently, using the isometry U, we choose S = L 2 (fl, p), 
so that the notation < designates an estimate in which the constant does not depend on u £ X, 
or on / £ L 2 (fl,p), but depends only on the mapping p. We will use these notations to provide 
estimates on the solutions u(-) of (13), meaning that the constants in the estimates do not depend 
on the solutions. 


Theorem 2 ([4]). In addition to the above assumptions, we assume that there exist T > 0 and 
Ct > 0 such that 

CW( 0 )||^ [ T \\BV*m\xdt, 

Jo 

for every solution of <j> f (t)+= 0 (observability inequality for the linear conservative equation). 

Then, for every uo £ X, there exists a unique solution u(-) £ C°(0, + 00 ; A')flC 1 (0, + 00 ; D(A)') 
of (13) such that u( 0) = uo- 2 Moreover, the energy of any solution satisfies 

E u (t) < Tmax( 7 i,£'„( 0 ))L Q,-i( 7z ^) ) ’ 

2 Here, the solution is understood in the weak sense (see [10, 18]), and D(A)' is the dual of D(A) with respect to 

the pivot space X. If uo E D(A), then u(-) E C°(0, -|-oo; D(A)) D C 1 (0, +oo; X). 
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for every time t ^ 0, with 71 ~ H-BH /72 arid, 72 — Ct/T(T 2 ||.B 1 / 2 || 4 + 1). If moreover 


limsup Atf(s) < 1, (14) 

sNjO 

then we have the simplified decay rate 

E u (t ) < Tmax( 7 l ,£ u (0)) (iA ) _1 (y) , 

for every time t ^ 0, for some positive constant 73 ~ 1 . 

Note the important fact that this result gives sharp decay rates (see Table 1 for examples). 


9(s) 

A h(s) 

decay of Eft) 

s/ ln p (l/s), p > 0 

lim sup A# (s) = 1 

e -t inp+1> / t i/(p+i) 

s p on [0, Sg], p > 1 

A h(s) = ^ < 1 

t -2/(p-l) 

e- 1 /^ 

lim Ah(s ) = 0 

s\,0 

1/ ln(t) 

s p ln 9 (l/s), p > 1, q > 0 


t -2/(p-l) ln -2 q /(p-i) ^ 

e -ln p (l/s) ; P > 1 

lim A h(s) = 0 

s\0 

e - 21 n 1 ^(i) 


Table 1: Examples 

Theorem 2 improves and generalizes to a wide class of equations the main result of [3] in which 
the authors dealt with locally damped wave equations. Examples of applications are given in [4], 
that we mention here without giving the precise framework, assumptions and comments: 

• Schrodinger equation with nonlinear damping (nonlinear absorption): 

idtu(t, x ) + A u(t, x ) + ib{x)u{t , x)p{x, | u(t, x)|) = 0 . 

• Wave equation with nonlinear damping: 

d u u(t , x) — A u(t, x) + b(x)p(x, d t u(t : x)) = 0. 

• Plate equation with nonlinear damping: 

d t tu(t, x) + A 2 u(t, x) + b(x)p(x, d t u(t, x)) = 0 . 

• Transport equations with nonlinear damping: 

d t u(t, x) + dW(v(x)u(t, x)) + b(x)p(x, u(t , x)) =0, ifT", 
with div(n) = 0 . 

• Dissipative equations with nonlocal terms: 

d t f + vS7 x f = p(f), 
with kernels p satisfying the sign assumption fp{f) ^ 0 . 
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Proof. It is interesting to quickly give the main steps of the proof. 

• Step 1. Comparison of the nonlinear equation with the linear damped model: 
Prove that the solutions of 


u\t) + Au(t) + BF(u(t )) = 0, 
z'(t) + Az(t) + Bz(t) = 0, z(0) = u(0), 


satisfy 

jj \\B 1/2 z(t)\\ 2 x dt < 2 £ (\\B l / 2 u(t)\\ 2 x + \\B 1/2 F(u(t))\\ 2 x ^j dt. 

• Step 2. Comparison of the linear damped equation with the conservative linear equation: 
Prove that the solutions of 


z'(t) + Az(t) + Bz{t) = 0, z(0)=u(0), 

+ Acj){t) = 0, 0(0) = u(0), 

satisfy 

[ \\B l/2 (t>{t)\\ 2 x dt ^k T f \\B 1/2 z{t)\\ 2 x dt, 

Jo Jo 

with k T = 8T 2 ||B 1 / 2 || 4 + 2. 

• Step 3. Following the optimal weight method introduced by F. Alabau (see, e.g., [2]), we 
set w{s) = L _1 (^j with /3 appropriately chosen. Nonlinear energy estimate: 

Prove that 

0 )) (\\B^ 2 u(t)\\ 2 x + \\B'/ 2 F(um 2 x) dt 

< T\\B\\H*(w{E40))) + {w{E 4> { 0)) + 1) [ {Bu(t ), F(u(t))) x dt. 

Jo 


• Step 4. End of the proof: 

Using the results of the three steps above, we have 

T\\B\\H*{w(E^0))) + (w(E^(0))^l) [ (Bu(t),F(u(t))) x dt 

Jo 

~ Jo w ( E ^^i} Bl/2u ^ 2x + ll- Bl/ 2 jF ( u W)ll^) dt (Step 3) 

> [ w(E (/> (0))\\B 1/2 z(t)\\ 2 x dt (Step 1) 

Jo 

>[ w{E c/) (0))\\B 1/2 4 l (t)\\ 2 x dt (Step 2) 

Jo 

> Cst w(E l j > (0))E(/,(()) (uniform observability inequality) 

from which we infer that 

E U (T) < E u { 0) (l - Pt L - 1 ££££) > 

and then the exponential decrease is finally established. 


□ 
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3.2 Space semi-discretizations 

In this section, we define a general space semi-discrete version of (13), with the objective of 
obtaining a theorem similar to Theorem 2, but in this semi-discrete setting, with estimates that 
are uniform with respect to the mesh parameter. As we are going to see, uniformity is not true 
in general, and in order to recover it we add an appropriate extra numerical viscosity term in the 
numerical scheme. 

3.2.1 Space semi-discretization setting 

We denote by Ax > 0 the space discretization parameter (typically, step size of the mesh), with 
0 < Ax < Ax 0 , for some fixed Ax 0 > 0. We follow [28, 29] for the setting. Let (.X’Ax)o<Aa;<Ax 0 
be a family of finite-dimensional vector spaces (X/\ x ~ with N{Ax) £ IN). 

We use the notations < and — as before, also meaning that the involved constants are uniform 
with respect to Ax. 

Let /? £ p{A) (resolvent of A). Following [18], we define Xi/ 2 = (/fid.v — A) _1 / 2 (X), endowed 
with the norm ||u||jc 1/2 = ||(/3idjc — A) l / 2 u\\x (for instance, if A 1 / 2 is well defined, then X\j 2 = 
Z^A 1 / 2 )), and we define X -i /2 = (dual with respect to X). 

The general semi-discretization setting is the following. We assume that, for every Ax £ 
(0,Ax o ), there exist linear mappings Pax '■ X_ l / 2 —> Xa x and Pax '■ Xa x —■* -^ 1/2 such that 
PaxPax = idx Aa ,i an d such that Pax = Pax- We assume that the scheme is convergent, that is, 
||(I — PaxPax)u\\x ~> 0 as Ax —> 0, for every u £ X. Here, we have implicitly used the canonical 
injections D(A) Xi/ 2 X X_x /2 (see [18]). 

For every Ax £ (0, Ax 0 ): 

• Xax is endowed with the Euclidean norm ||uax||ax = \\PaxUax\\x, for uax 6 Xax- The 
corresponding scalar product is denoted by (v)ai 3 

• We set 4 _ _ 

Aax — PaxAPaxi P Ax = PAxB PAx • 

Since Pax = Pa x > Aax is skew-symmetric and Ba x is symmetric nonnegative 

• We define F Ax ■ X Ax -» X Ax by 

VuAx £ Aa x FAx(FAx) = PA X F [PAx^ Ax) 

Note that Bax is uniformly bounded with respect to Ax, and Fax is Lipschitz continuous on 
bounded subsets of Xa x , uniformly with respect to Ax. 

Now, a priori we consider the space semi-discrete approximation of (13) given by 

^Ax(^) T dAx^Ai(i) T B Ax Pax {FAx (t)) ~ dj (15) 

and we wonder whether we are able or not to ensure a uniform decay rate of solutions of (15). 

It happens that the answer is NO in general. Let us give hereafter an example in the linear 
case (F = idv - )- 

3 Note that \\Pax\\l(x &x ,x) = 1 and that, by the Uniform Boundedness Principle, \\Pax\\l(x,x A x ) iS 1- 

4 We implicitly use the canonical extension of the operator A : X\ /2 -A X_i/ 2 - 
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3.2.2 Finite-difference space semi-discretization of the ID damped wave equation 

Let us consider the ID damped Dirichlet wave equation 


ytt - Vxx + ay t = 0 , 0 < x < 1 , t > 0 , 

2 / 0 , 0 ) = 2 / 0 , !) = 0 , 


(16) 


with a ^ 0 measurable and bounded, satisfying a(x) )a>0onwC (0,1) of measure 0 < |w| < 1. 
This equation fits in our framework with 


The energy along any solution is 


A = 


0 id 
-A 0 


B = a id, F = id. 


and we have 


1 f 1 

E(t) = -J (y t (t,x) 2 +y x (t,x) 2 )dx, 
E'(t) = — / a(x)y t (t,x) 2 dx < 0. 

J o 


It is known that Eft) has an exponential decrease. Indeed, by Theorem 1, the exponential decrease 
is equivalent to the observability inequality 


f 

! 0 Juj 


0tO,x ) 2 dx dt ^ C f (0 t (O, x) 2 + 0(0, x) 2 ) dx 

• Jo 

for every solution 0 of the conservative equation 


(17) 


<t>tt ~ <t>xx = o, 00 , 0 ) = <j>(t, 1 ) = 0 , (18) 

and we have the following result. 

Proposition 1. The observability inequality (17) holds, for every T ^ 2. 

Proof. The proof is elementary and can be done by a simple spectral expansion. Any solution 0 
of (18) can be expanded as 


00 ,®) = 


OO 

E 

k—l 


a/c cos(/c7rt) + E- sin( knt) ) sin(/c7rx), 

K7T 1 


with 


Then, 


/o 

Now, for every T ^ 2, we have 

rT 


0 ( 0 , x) = E a k sin(fc 7 rx), and 0 t (O, x) = b k sin(fc7rx). 
fc=1 fc=l 

A i °° 

(0t(O, a ;) 2 + 0(0, x) 2 ) dx = -^(a 2 +bj:). 


k—l 


' 0 J uj 


(ft{t,x) 2 dx dt ^ / / (j> t {t,x) 2 dxdt 

J 0 ui 
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and 


,•2 [ oo f 2 

/ / (f> t (t,x) 2 dxdt = ^ / (aj cos(jnt) + bj sm(jnt))(cik cos(fc 7 rf) + &*, sin(fc 7 rf)) dt 

J 0 J oj j fc = i ^ 0 

x / sin(j 7 ra;) s\n(kiTx) dx 
J (jj 

OO „ 

= 5Z( a i + ^j) / sin 2 O' 7 ™) dx 

J OJ 


At this step, we could simply conclude by using the fact that the sequence of functions sin 2 (j7rx) 
converges weakly to 1/2, and then infer an estimate from below. But it is interesting to note that 
we can be more precise by showing the following simple result, valid for any measurable subset u> 
of [0,1]. 

Lemma 2. For every measurable subset w C [0,1], there holds 


,. [ . n,. . , U - sin U 

\/j £ Iht / sm 2 (j7rx)dx > ----, 

J OJ " 

where |w| is the Lebesgue measure of u>. 

With this lemma, the observability property follows and we are done. It remains to prove the 
lemma. 

For a fixed integer j , let us consider the problem of minimizing the functional 



sin 2 (jiTx) dx 


over the set of all measurable subsets oj' C [0,1] of Lebesgue measure \uj'\ = |w|. Clearly there 
exists a unique (up to zero measure subsets) optimal set, characterized as a level set of the function 
x i ^ sin 2 (jirx), which is the set 


U! 




and there holds 


r ^ r M/ 2 J r\u\r* ^ 

/ shr(jTrx) dx = 2 j / sin 2 jxdx = 2 / sin 2 udu = -(|w| — sin(|w|)). 

A‘ nf Jo Jo 2 

Since the quantity is independent of j , the lemma follows, and thus the proposition as well. 


M/2 


□ 


Hence, at this step, we know that the energy is exponentially decreasing for the continuous 
damped wave equation. 

Let us now consider the usual finite-difference space semi-discretization model: 


Vi- 


Vj+1 - 2 Vj + Vj-1 
(Aa:) 2 


+ ajyj = 0, 


2/o (i) = VN+i{t) = 0. 


= 1, 


, AT, 
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The energy of any solution is given by 




and we have 


N 


e 'a x (t) = ~ Ax J2 a ^y'j^ 2 ^ °- 

3 =1 

We claim that E^ x (t) does not decrease uniformly exponentially. 
Indeed, otherwise there would hold 


N r T 


E Ax {0) < CAx Y, 

3=1 ' 

for every solution of the conservative system 

2 <t>j +1 — + cj)j -1 

** {Ax) 2 

<l>o(t) = </>N+l(t) = 0, 


a j{^{t)) 2 dt 


= 0, 




(19) 


for some time T > 0 and some uniform C > 0. But this is known to be wrong. To be more precise, 
it is true that, for every fixed value of Ax, there exists a constant Cax , depending on Ax, such 
that any solution of (19) satisfies 

N r T 

E Ax ( 0 ) < C Ax Ax Y I aj{^j{t)) 2 dt. 

3=1 

This is indeed the usual observability inequality in finite dimension, but with a constant C Ax 
depending on Ax (by the way, this observability inequality holds true for any T > 0, but we are 
of course interested in taking T ^ 2 in a hope to be able to pass to the limit). Now, the precise 
negative result is the following. 

Proposition 2. [33] We have 

Cax —> +oo. 

Ate—>•() 

In other words, the observability constant is not uniform and observability is lost when we pass 
to the limit. The proof of this proposition combines the following facts: 

• Using gaussian beams it can be shown that along every bicharacteristic ray there exists a 
solution of the wave equation whose energy is localized along this ray. 

• The velocity of highfrequency wave packets for the discrete model tends to 0 as Ax tends to 

0. 


Then, for every T > 0, for Ax > 0 small enough there exist initial data whose corresponding 
solution is concentrated along a ray that does not reach the observed region oj within time T. The 
result follows. 


Remark 4. This proof of nonuniformity is not easy in the case of an internal observation, and 
relies on microlocal considerations. It can be noticed that nonuniformity is much easier to prove, 
in an elementary way, in the case of a boundary observation. We reproduce here the main steps 
of [24] with easy-to-do computations: 
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1. The eigenelements of the (finite-difference) matrix Aai are 


( \ 


x,j — 


2. If Ax —> 0 then 


sin(j7rfcAx) 


V 


/ 


4 . 2 jirAx 

(Ax) 2 


A A x,i = -7*^2 Sin 2 j = 1 , • • •, AT. 


\j |AAx,n| ^ | Aax,ti— 1 1 ^ Ax. 

3. The vector-valued function of time 

^A*(f) = , (e lt ^^*- N 'e Ax , N - 

\/|Aax,]v| v ' 

is solution of </>a x(t) = Aai^AiW, and if Ax tends to 0 then 

N r u j+1 ( 0 )-^( 0)| 2 


the term Ax 

3=0 

the term Ax 


/ <j>N(t) \ 


Ax 

i 


1^(0) | has the order of 1, 


/o \ Ax J 

4. Therefore Ca x —> +oo as Ax —> 0. 


dt has the order of Ax. 


Remedy for this example. A remedy to this lack of uniformity has been proposed in [40], 
which consists of adding in the numerical scheme a suitable extra numerical viscosity term, as 
follows: 


Vi 


Vj +1 - 2 Vj + Vj _i . 2 ( Vj +1 ~ 2 Vj + Vi -1 


(Ax) 2 
2/oW = 2 /jv+i(*) = 0. 


+ - ( A ^)" 


\ (Ax) 2 


= 0, j = 


With matrix notations, this can be written as 

2/Ax A Ax2/Ax T ^AxZ/Ax (Ax) A Ax2/Ax = 0. 
Now the energy decreases according to 


AT 


N 


Eax = -Ax ^ dji/j - (Ax) 3 ^2 


a=i 


3=0 


( 2/j+i - Vj 

\ Ax 


and we have the following result. 

Proposition 3 ([40]). There exist C > 0 and 5 > 0 such that, for every Ax, 

EAx(t) ^ C exp(— <5/)£ i ax(0). 
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In other words, adding an appropriate viscosity term allows to recover a uniform exponential 
decrease. 

An intuitive explanation is that the high frequency eigenvalues are of the order of ■ Indeed, 


X 2 

A Ax,j 


(Ax)‘ 


■ Sill 


j tt Ax 


= 1 , 


,N. 


Then for highfrequencies the viscosity term (Ax) 2 A Ax y Ax has the order of ji^, and hence it 
behaves as a damping on highfrequencies. The main result of [40] shows that it suffices to recover 
an uniform exponential decay. 

Their proof mainly consists of proving the uniform observability inequality 


N r T 

E Ax ( 0) i CAx J2 


2 I I ‘Ai+l 


a. j(<Pj(t)) +(Ax) 


J=i J0 


Ax 


dt , 


for every solution of 

<Aj+i + 2 <f>j A <f>j- 1 4>j+i — 2<t>j A <t>j- 1 

4 (Ar) 2 

<Ao (t) = <pN+i(t) = 0, 


= 0, j = 1,...,N, 


and this follows from a meticulous application of a discrete version of the multiplier method. 

Note that they also prove the result in a 2D square, and that this was generalized to any regular 
2D domain in [35]. 


3.2.3 General results for linear damped equations, with a viscosity operator 

The idea of adding a viscosity in the numerical scheme, used in [40] in order to recover a uniform 
exponential decay rate for the finite-difference space semi-discretization of the linear damped wave 
equation, can be generalized as follows. 

We consider (13) in the linear case F = id_\', that is, 

u'(t) + Au(t ) + Bu(t) = 0, 

with A : D(A) —> X skew-adjoint and B = B* ^ 0 bounded. The energy along any solution is 

and we have 

E(t) = -WB^um 2 ,. 

We consider the space semi-discrete model with viscosity: 

A A/\ x u/\ x A + (Ax) = 0} (20) 


where Vax : X& x —> X& x is a positive selfadjoint operator, called the viscosity operator. 

In the previous example (ID damped wave equation), the viscosity term was —(Ax) 2 A& x y' Ax , 
hence 


a = 2, 


V A * 


(° ° ) 
Vo -a A x .) 


Note that it is not positive definite, however it was proved to be sufficient to infer uniformity. 
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In general, as we are going to see, a typical example of a viscosity operator is Vax = y 

plus, possibly, eidx Ax (to make is positive). 

In any case, the energy along a solution of (20), defined by 

(0 = 7^||uAx(t)|| Ail 

has a derivative in time given by 

E UAx (t) = — 11 Ax u Ax (t) 11 Ax ^ ( A x) a \\(V Ax ) 1 / 2 u Ax (t)\\ 2 Ax < 0 . 

Now, we have the following result. 

Theorem 3 ([21]). Assume that the uniform observability inequality 

cu Ax m 2 < £ + (Axywv^^m 2 ) dt 

holds for every solution of the conservative system 

pAx T A/^ x cf)/\ x = 0. 


Then the uniform exponential decrease holds, i.e., there exist C > 0 and S > 0 such that 


/Sx ( t ) ^ C exp(—<5t)£' UAx (0), 


for every Ax, for every solution mai of (20). 

This theorem says that, in general, adding a viscosity term in the numerical scheme helps to 
recover a uniform exponential decay, provided a uniform observability inequality holds true for 
the corresponding conservative equation. Of course, given a specific equation, the difficulty is 
to establish the uniform observability inequality. This is a difficult issue, investigated in some 
particular cases (see [48] for a survey). 

Let us mention several other related results. In [19], such uniform properties are established for 
second-order equations, with a self-adjoint positive operator, semi-discretized in space by means of 
finite elements on general meshes, not necessarily regular. The results have been generalized in [34]. 
In [1], one can find similar results under appropriate spectral gap conditions (in this reference, we 
can also find results on polynomial decay, see further). In [36, 37], viscosity operators are used for 
the plate equation, and for general classes of second-order evolution equations under appropriate 
spectral gap assumptions, combined with the uniform Huang-Priiss conditions derived in [32]. 

Remark 5. In order to recover uniformity, there exist other ways than using viscosities (see [48]: 

• Use an appropriate numerical scheme (such as mixed finite elements, multigrids). 

For instance, considering again the ID damped wave equation (16), the mixed finite element 
space semi-discretization model is 


Vj+i + 2 Vj + Vj-i Vj+i ~ 2 Vj + Vj -1 


(Az) 2 


yo(t) = yN+i{t) = o, 


J 4 1/2 fa-i + Vj) + Qj ^ 1/2 (Vj + Vj+ 1) = 0, 


The energy of any solution is 

(( fe+ l(t) 2 +feW ) a + (AiAz jATQ , 
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and one has 


N—l 

EAx{t) = - Ax '^2 a j + 1/2 
3=0 


Vj+i(t) +yj(t) 


€ 0. 


It is then proved in [20] that Ea x {4) decreases exponentially, uniformly with respect to Ax. 
The proof consists of proving that the uniform observability inequality: 


N r T 


-E’ax(O) ^ CAx 'y ] f a,j_ |_i /2 

3=1 J ° 


^j+iif) + Mt) 


dt , 


for any solution of 


2 4 * 3+1 - 2 < A? + <t>j-i 

% (Ax) 2 

4>o(t) = 4 >n+i (t) = o. 


= 0 , 


j = 1,...,N, 


The proof of this uniform observability inequality is done thanks to a careful spectral analysis, 
using Ingham’s inequality. 

The main difference with the finite-difference numerical scheme is that, here, the discrete 
eigenvalues are given by 


X 2 

^Aaij 


(Ax) 


■ tan' 


jnh 


= 1 , 




and then their behavior at highfrequencies ( j close to N) is very different from the case of 
finite differences where the tan function was replaced with the sin function: the derivative of 
the discrete symbol is now bounded below and hence there are not anymore any wavepackets 
with vanishing speed (see [48]). 


• Filtering highfrequencies. 

In the general framework 

u + Au + Bu = 0, 

with A : D(A) -+ X skew-adjoint and B = B* ^ 0 bounded , considering the semidiscrete 
model 

ilAx T A Ax’ll Ax + B Ax’ll Ax — 0, 

the idea (well surveyed in [48]) consists of filtering out the highfrequencies, so that the uniform 
observability inequality 

C'||^ £c (0)|| 2 < f \\B l H<fAx{t )\\ 2 dt 
Jo 


holds for every solution of the conservative system cj>Ax — Aax4*Ax with 


4*Ax (0) £ Cs/(Ax)”(Aax) (spectrally filtered initial data), 


with 


Cs^Aax) — SpanjeAzj | |Aax,j| ^ a}- 


Here, the vectors eAx,j are eigenvectors of Aa x associated with the eigenvalues XaxJ • Then 
the uniform exponential decrease holds. 


Actually, as noticed in [21], we have equivalence of: 
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— Every solution of cf>Ax = Aai^Ai satisfies 

ci||0A*(o)|| 2 < £ (wB^Axim 2 + (AxywvgMm 2 ) dt. 

- Every solution of <pAx = Aax^Ax with 0 a® (0) G C^As^C^a®) satisfies 

c 2 UAxm 2 < 

Jo 

In conclusion, there are several ways to recover uniform properties for space semi-discretizations: 

1. Add a viscosity term (sometimes called Tychonoff regularization): its role being to damp out 
the spurious highfrequencies. 

2. Use an adapted method (like mixed finite elements, or multi-grid methods), directly leading 
to an uniform observability inequality. 

3. Filter out highfrequencies, by using spectral projections. 

In our case, we are going to focus on the use of appropriate viscosity terms, in order to deal more 
easily with nonlinear terms. 

3.2.4 Viscosity in general semilinear equations 

We now come back to the study of the general semilinear equation (13), and according to the 
previous analysis, we consider the space semi-discrete approximation of (13) given by 

u' Ax (t) + A A xUAx{t) + B Ax FAx{uAx{t)) + (Ax) a VAxUAx{t) = 0. (21) 

As before, the additional term (Ax) a VAxUAx{t), with a > 0 , is a numerical viscosity term whose 
role is crucial in order to establish decay estimates that are uniform with respect to Ax. We only 
assume that Vax ■ Xax —> Xa x (viscosity operator) is a positive selfadjoint operator. 

Our objective is to be able to guarantee a uniform decay, with the decay rate obtained in 
Theorem 2. 

The energy of a solution uax of (21) is 

Fu &x (t) = — ||UAa:(i)|| An 

and we have, as long as the solution is well defined, 

^« A ,W = ~{ u Ax{t), BAxFAx{uAx{t))) Ax ~ (Ax) ||(Vax) ^ u Ax(t)\\ a x - 
We assume that 

(uAx,Ba x Fax(uAx))ax + {.Axy\\(VAx) 1/2 UAx\\ 2 A x > 0 , 

for every uax G Xax- Hence the energy E u&x (t) is nonincreasing. 

For every / £ T 2 (H, /x), we set 

pAx(f) = U^pAxPAxUptf) = U~^ PAx FAx F(U f). 

The mapping pAx is the mapping p filtered by the “sampling operator” 

U-'PaxPaxU = (PaxU)*PaxU. 
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By assumption, the latter operator converges pointwise to the identity as Ax —> 0, and in many 
numerical schemes it corresponds to take sampled values of a given function /. 

We have /5 ax( 0) = 0. Setting f a x = U~ 1 PaxUa x , we assume that 

fAxpAx(fAx) ^ 0 , 

and that 

Cl g{\fAx(x)\) < \pAx{fAx)(x)\ ^ c 2 g~ 1 (\fAx(x)\) for a.e. x <E Vt such that |/ax(^)| < 1, 

Cl |/ax(^)| < \pAx{fAx)(x)\ ^ C 2 |/Ax(a;)| for a.e. i£l! such that |/ax(z)| > 1, 

for every uax € ^Ai, for every Ax € (0 ,Axq). 

Note that these additional assumptions are valid for many classical numerical schemes, such 
as finite differences, finite elements, and in more general, for any method based on Lagrange 
interpolation, in which inequalities or sign conditions are preserved under sampling. But for 
instance this assumption may fail for spectral methods (global polynomial approximation) in which 
sign conditions may not be preserved at the nodes of the scheme. 

Theorem 4 ([4]). In addition to the above assumptions, we assume that there exist T > 0, a > 0 
and Ct > 0 such that 

CtE^M € £ (lli^AxCOIlL + (AoOIv^axWIIax) dt , (22) 

for every solution of <p' Ax (t) + AAx4>Ax(t) = 0 (uniform observability inequality with viscosity for 
the space semi-discretized linear conservative equation). 

Then, the solutions of (21), with values in Xa x , ere well defined on [0, +oo), and the energy 
of any solution satisfies 

E u^{t) <Tmax(7i,£; UAa ,(0))L £_££ , 

for every t ^ 0, with 71 ~ H-BH /72 end 72 ~ Ct/T(T 2 ||13 1/,2 || 4 + 1 ). Moreover, under (14), we 
have the simplified decay rate 

Eu^it) < Tmax( 7 l ,£ UAi ( 0 )) [H')~ l (y) , 

for every t ^ 0 , for some positive constant 73 ~ 1 . 

The main assumption above is the uniform observability inequality (22), which is not easy 
to obtain in general, as already discussed. We stress again that, without the viscosity term, 
this uniform observability inequality does not hold true in general. As said previously, a typical 
example of a viscosity operator, for which uniform observability results do exist in the literature, 
is Vax = A^ x Aa X ‘ 

3.3 Time semi-discretizations 

In this section, we show how to design a time semi-discrete numerical scheme for (13), while keeping 
uniform decay estimates. 

As a first remark, let us note that a naive explicit discretization 

/c - j-1 /c 

————-—b Au k + BF(u k ) = 0, tk = kAt, 

is never suitable since 
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• it is instable and does never satisfy the basic CFL condition (note also that we are going to 
consider full discretizations in the next subsection), 


• it does not preserve the energy of the conservative part. 

Hence an implicit discretization is more appropriate. 
Following [21], we choose an implicit mid-point rule: 


u k+1 - u k 
At 


A 


u k + u k+1 


BF 


u k + u k+1 


= 0 , 


tk = kAt. 


Defining the energy by 


we have 

E u k +1 — E u k = —At 



+ u k+1 ^ 


,F 




X 


The question is then to determine whether we have or not a uniform decay. 

As before, the answer is NO in general, and in order to see that, we are first going to focus on 
the linear case. 


3.3.1 Linear case 

Assuming that F = id_v, we consider the equation 

v! + Au + Bu = 0, 


with the time semi-discrete model 


u k+1 - u k 
A t 


A 


We have 


E,,k = -llw' 




u k + u k+1 


E~.k+1 — 


B 


u k + u k+1 


= 0 . 


B 1/2 


u k + u k+1 


< 0 , 




2" At 

and we investigate the question of knowing whether we have a uniform exponential decrease 

E u k ^ CE U o exp(-SkAt). 


A negative answer is given in [46] for a linear damped wave equation. We do not provide the 
details, that are similar to what has already been said. 

Therefore, as before, filtering, or adding an appropriate viscosity term, is required in a hope to 
recover uniformity. 

In [21], the authors propose to consider the following general numerical scheme: 


~k+i _ u k 

At 

u k+1 - u k+1 

At 


A 


u k + u k+1 


B 


u k + ~k+1 


= 0 , 


= V A tu k+1 , 
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where Va t ■ X —> X ia s positive selfadjoint operator (viscosity operator). Written in an expansive 
way, we have 


u k+1 _ u k 


At 


A 


u k + u k+1 


B 


,fe+i 


At 


At 


+ —SV At u fe+1 + V At u k+1 + — AV At u k+1 = 0. 


Considering the energy E u k = 2 ||ti fc ||x, we have 


E„ 


— E lt k 


At 


B 1 / 2 


'lk -\-1 


- \\(V At ) 1/2 u k+1 \\ 2 x - 


At 


||V At M fe+1 ||^ <0. 




Now, the main (remarkable) result of [21] is the following. 

Theorem 5 ([21]). 1. We assume the uniform observability inequality with viscosity for the 

time semi-discretized linear conservative equation with viscosity: 


N -1 


CtUYx < a tJ2 


k—0 


B 


1/2 




2 

X 


N -1 TV—1 

+ AtJ2 ll(VA t ) 1/2 ^ +1 ||i + (Af) 2 ^ ||W +1 ||i 

fe =0 fc =0 


for all solutions of 




(j) k+l - i 

A T 

cj) k+1 - i 

At 


, fc 

- + A 

,k+1 



= 0 


Then the uniform exponential decrease holds, i.e., 

E u k ^ CE U o exp(— 5kAt), 

for some positive constants C and 6 not depending on the solution. 
2. The previous uniform observability inequality holds true if 

V At = -(At) 2 A 2 = (At) 2 A* A. 


Note that this result provides a choice viscosity operator that always works, in the sense that 
the uniform observability inequality is systematically valid with this viscosity. Other choices are 
possible, such as 

V A * = — (id.v - (At) 2 A 2 )- 1 (At) 2 A 2 . 

The proof is done by resolvent estimates (Hautus test) and by carefully treating lower and higher 
frequencies. 
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3.3.2 Semilinear case 


Coming back to the general semilinear case, we consider the implicit midpoint time discretization 
of (13) given by 


u k+1 - u k 
At + 

~k+l _ u k +1 

At 

, u° = u( 0 ). 


u k + u k+1 


+ BF 


u k + u k+1 


= 0 , 


= V At u k+1 , 


(23) 


Written in an expansive way, (23) gives 


u k+1 - u k 
At 


+ A 


u k + u k+1 


BF 


+ (id A - + A tV At )u k+1 


V At u k+l + ^ AV At u k+1 = 0. 


For the energy E u k = 2 W™ I I jo we have 


E u k + 1 — E u k = -At (\B 


u k + u k+1 


,F 


u k + u k+1 


X 


- At \\{V At f/ 2 u k+1 f x - ^\\V At u k+1 \\ 2 x < 0, 

for every integer k. 

Theorem 6 ([4]). In addition to the above assumptions, we assume that there exist T > 0 and 
Ct > 0 such that, setting N = [T /At] (integer part), we have 


N -1 




k—0 


B !/2 ( ^ + / +1 ' 


X 


N—l 


N—l 


+ At ^ ||(V At ) 1/2 ^ +1 ||^ + (At) 2 £ \\V At <t> k+1 fx, 


k -0 


k—0 


for every solution of 


■kk+l 


ik , Ak+1 


At 

A+l _ +k+l 


= 0 , 


At 


= V At cf k+ k 


(uniform observability inequality with viscosity for the time semi-discretized linear conservative 
equation with viscosity). 

Then, the solutions of (23) are well defined on [0, +oo) and, the energy of any solution satisfies 


E u k "A T max(yi, E u o) L 


1 


ip~ 1 ('Y2kAt) 


for every integer k, with 72 ~ CV/T(1 +e 2T ^ B ^ max(l, T||B||)) and 71 ~ 11T?11 / 72 . Moreover, under 
(14), we have the simplified decay rate 

E uk < T max( 7 i , E u o ) (iT)” 1 ' 73 

for every integer k, for some positive constant 73 ~ 1 . 


kAt J ’ 
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3.4 Full discretizations 


Results for full discretization schemes can be obtained in an automatic way, from the previous 
time discretization and space discretization results: indeed, following [ 21 ], we note that the results 
for time semi-discrete approximation schemes are actually valid for a class of abstract systems 
depending on a parameter, uniformly with respect to this parameter that is typically the space 
mesh parameter Ax. Then, using the results obtained for space semi-discretizations, we infer the 
desired uniform properties for fully discrete schemes. In other words, we first discretize in space 
and then in time. We refer to [4] for the detailed definition of that abstract class. 

Theorem 7 ([4]). Under the assumptions of Theorems 2, 4 and 6, the solutions of 


u k +} — u‘ 


At 

fc+1 ,, k +1 


A r / U 

— + a Ax 


~k +1 ' 


Ai w “Aj 


T E Ax F Ax 


A I y k + l 
l Ax + U Ax 


rfc+1' 


Va, 


‘Ai w “Ai 


= 0 , 


u Ax - u 


At 


Ax fc+1 

= Va tu£ x , 


are well defined for every integer k, for every initial condition u Ax £ X Ax , and for every Ax £ 
( 0 , Axo), and the energy of any solution satisfies 

^aJIL = E <x < , 

for every integer k, with 72 ~ Ct/T( 1+ max(l, T||R||)) and 71 ~ ||R||/ 72 . Moreover, under 

(14), we have the simplified decay rate 


kAt J 

for every integer k, for some positive constant 73 ~ 1 . 

As an example, we consider the nonlinear damped wave equation 


duuit , x) — Au(t, x) + b(x)p{x , dtu{t , x)) = 0 , 


with appropriate assumptions on p , as discussed previously. We first semi-discretize in space with 
finite differences with the viscosity operator V Ax = ~(Ax) 2 A Ax , obtaining 

u Ax,a(f) A Ax u Ao0}Cr (t) + b Ax (x a ) p(x a , n Axa (t)) (Ax) A Ax u Ax (t) = 0, 


where a is an index for the mesh. The solutions of that system have the uniform energy decay 
rate L(\/if~ l (t)) (up to some constants). 

We next discretize in time, obtaining 


-•fc+i 


- v 


fc+i 
A Xi<T 


At 

n .k 



At 

~ k +l 

..fc+i 

a Ax,o 

U Ax,a 


At 


„k+l 

^A x,cr 

^A x,cr 


U Ax,a + U Ax]a lU \„( „ V Ax,cr + ^Ax.a 

^Ax,(T 0 ^Arc,cr \%cr )P I %cr: n 


- (Ax) 2 Aa x,a 




= 0, 


^=-(At) 2 ALxt: 


= -(AtfA 2 v k + x \, 


At 
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and according to Theorem 7, the solutions of that system have the uniform energy decay rate 
L{\/xf~ 1 {f)) (up to some constants). 

3.5 Conclusion and open problems 

3.5.1 (Un)-Boundedness of B 

We have assumed B bounded: this involves the case of local or nonlocal internal dampings, but 
this does not cover, for instance, the case of boundary dampings. 

For example, let us consider the linear ID wave equation with boundary damping 

dttu - d xx u = 0, t G (0, +oo), x G (0,1), 

u(t, 0) = 0, d x u(t, 1) = — ad t u(t, 1), fG(0,+oo), 

with a > 0. The energy 

E(t) = ^ [ (\d t u(t,x)\ 2 + \d x u(t,x)\ 2 ) dx 

^ J o 

decays exponentially. It is proved in [41] that the solutions of 

U^ x (Aj?) ^As^Aj: = 

UAx{t, 0) = 0, D Ax UAx(t, 1) = - ad t u(t , 1), 

(regular finite-difference space semi-discrete model with viscosity) with 

1 ... 0 \ 

•• 1 
.. 0 -l) 

(forward Unite-difference operator), have a uniform exponential decay. But this fails without the 
viscosity term. 

This example is not covered by our results. 

We also mention [6] where a sufficient condition (based on energy considerations and not on 
viscosities) ensuring uniformity is provided and is applied to ID mixed finite elements and to 
polynomial Galerkin approximations in hypercubes. The extension to the semilinear setting is 
open. 

3.5.2 More general nonlinear models 

If the nonlinearity F involves an unbounded operator, like for instance 

u'(t) + Au(t) + BF(u(t),Vu(t)) = 0, 
then the situation is widely open. 

For instance, one can think of investigating semilinear wave equations with strong damping 
d t tu — A u — aAd t u + bd t u + g* A u + f(u) = 0, 

with / not too much superlinear. There are many possible variants, with boundary damping, with 
delay, with nonlocal terms (convolution), etc. Many results do exist in the continuous setting, but 
the investigation of whether the uniform decay is preserved under discretization or not remains to 
be done. 


/- 1 

0 

u 
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3.5.3 Energy with nonlinear terms 


Another case which is not covered by our results is the case of semilinear wave equations with 
locally distributed damping 

d u u — A u + a{x)d t u + f(u ) = 0, 

with Dirichlet boundary conditions, a a nonnegative bounded function which is positive on an 
open subset w of fi, and / a function of class C 1 such that /(0) = 0, sf(s) > 0 for every s £ 1R 
(defocusing case), |/'(s)| < C|s| p_1 with p < n/(n — 2) (energy subcritical). It is proved in [47] 
(see also [17, 25]) that, under geometric conditions on u>, the energy involving a nonlinear term 

[ (( d t u ) 2 + || Vu || 2 + F(u)) dx, 

J n 

with F(s) = Jq f, decays exponentially in time along any solution. The investigation of the uniform 
decay for approximation schemes seems to be open. In particular, it would be interesting to know 
whether microlocal arguments withstand the discretization procedure. In particular we raise the 
following informal question: 

Do microlocalization and discretization commute? 

3.5.4 Uniform polynomial energy decay without observability 

Let us assume that F = idx in (13) (linear case). As proved in Theorem 1, observability for the 
conservative linear equation is equivalent to the exponential decay for the linear damped equation. 

If observability fails for the conservative equation, then we do not have an exponential decay, 
however we can still hope to have, for instance, a polynomial decay. This is the case for some 
weakly damped wave equations when GCC fails (see [9]). 

In such cases, it is interesting to investigate the question of proving a uniform polynomial 
decay rate for space and/or time semi-discrete approximations. In [1], such results are stated for 
second-order linear equations, with appropriate viscosity terms, and under adequate spectral gap 
conditions. The extension to a more general framework (weaker assumptions, full discretizations), 
and to semilinear equations, seems to be open. 
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